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Abstract 

A two-dimensional computational fluid dynamics model is used to predict the oscillatory 
flow through a tapered cylindrical tube section (jet pump) placed in a larger outer tube. 
Due to the shape of the jet pump, there will exist an asymmetry in the hydrodynamic end 
effects which will cause a time-averaged pressure drop to occur that can be used to cancel 
Gedeon streaming in a closed-loop thermoacoustic device. The performance of two jet pump 
geometries with different taper angles is investigated. A specific time-domain impedance 
boundary condition is implemented in order to simulate traveling acoustic wave conditions. 
It is shown that by scaling the acoustic displacement amplitude to the jet pump dimensions, 
similar minor losses are observed independent of the jet pump geometry. Four different flow 
regimes are distinguished and the observed flow phenomena are related to the jet pump 
performance. The simulated jet pump performance is compared to an existing quasi-steady 
approximation which is shown to only be valid for small displacement amplitudes compared 
to the jet pump length. 


PACS numbers: 43.35.Ud, 43.25.Nm, 43.20.Mv, 47.32. Ff 
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FIG. 1: Jet pump with parameters that dehne the geometry (not to scale). Bottom dashed 
line indicates centerline, top solid line indicates tube wall. 


I. INTRODUCTION 

A jet pump is a crucial part of most closed-loop thermoacoustic devices.^ In such devices, a 
time-averaged mass flux known as Gedeon streaming can exist. ^ This time-averaged mass flux 
results in convective heat transport that can severely degrade the efficiency of thermoacoustic 
devices.^ To suppress Gedeon streaming, a jet pump can be used. Backhaus and Swift have 
shown that by correctly shaping a jet pump it is possible to take advantage of asymmetric 
hydrodynamic end effects to impose a pressure drop across the jet pump.® A typical jet pump 
consists of a narrowed, tapered tube section as shown schematically in Fig. [T] By balancing 
the pressure drop across the jet pump with that which exists across the regenerator, it is 
possible to produce a net zero time-averaged mass flux in the thermoacoustic device. 

Despite the proven effectiveness of jet pumps, there is a lack of understanding with 
respect to the exact fluid dynamics that lead to the observed pressure drop. Gurrent criteria 
for the design of a jet pump assume that the flow at any point in time has little “memory” 
of its past history — which is often referred to as the Iguchi-hypothesis.® This allows the 
acoustic behavior to be based on a quasi-steady approximation using minor loss coefficients 
reported for steady pipe flow.® 


A. Quasi-steady theory 


The pressure drop generated by an abrupt pipe transition in steady flow can be calculated 
using 


^Pml 



( 1 ) 


where K is the minor loss coefficient which depends on geometry and flow direction, p is 
the fluid density and u is the fluid velocity. For an abrupt expansion, K = Kg^p and can be 
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estimated using the Borda-Carnot equation.l^ In the case of a uniform flow 

Ke.p = (l - . (2) 

where Ag is the smaller cross-sectional area before the expansion and Aq is the cross-sectional 
area right after the expansion. Note that these values assume a uniform velocity prohle. Non- 
uniform velocity prohles will result in larger minor loss coefficient values for expansion.^ 
For a contraction, the steady flow minor loss coefficient is dependent upon the dimen¬ 
sionless curvature of the transition, Rc/D, where Rc is the radius of curvature and D is the 
diameter of the opening.l^For a sharp contraction, Kcon = 0.5 but this reduces to Kcon = 0.04 
for Rc/D > 0.15. 

Under the assumption that the Iguchi-hypothesis is applicable and that the minor loss 
coefficients have the same values in oscillatory flow as they do in steady flow, a quasi-steady 
model has been formulated by Backhaus and Swift to calculate the time-averaged pressure 
drop across a jet pumpP 


^P2,JP 
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(3) 


where is the velocity amplitude at the small exit of the jet pump. The subscripts “s” 

and “6” indicate the small and big opening of the jet pump, respectively. 

Although this time-averaged pressure drop can be exploited to cancel Gedeon streaming 
and improve the efficiency of a looped thermoacoustic device, this approach is not without 
penalty. Adding a jet pump results in additional dissipation of acoustic power. Under the 
same previous assumptions, the time-averaged acoustic power dissipation across a jet pump 
isSl 
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Qualitative evidence exists which supports the current analysis, but quantitative agree¬ 
ment between the theory and experiments remains poor.EEl While the accuracy of this ap¬ 
proach is yet unknown, it is assumed valid for large displacement amplitudes in relation to 
the jet pump dimensions.® Moreover, when using minor loss coefficients for steady expan¬ 
sion and contraction, the effect of the jet pump taper angle or the jet pump length is not 
included in the current theory while it is observed to have an important effect on the jet 
pump pressure drop.®' 
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B. Literature review 

Previous studies related to jet pumps for thermoacoustic applications include mainly ex¬ 
perimental or applied work; only a few computational studies have been published to date. 
Petculescu and When measured the pressure drop for a series of jet pump geometries in a 
standing wave experimental apparatus.® They then derived minor loss coefficients based on 
the measured pressure and the velocity in the jet pump waist which was estimated using 
an acoustic network model. A difference between the measured and theoretical minor loss 
coefficients is reported, especially for the diverging flow direction. Nevertheless, for the in¬ 
vestigated geometries - up to a taper angle of 10° - good agreement between the performed 
steady flow and oscillating experiments is obtained. An increase in the taper angle is shown 
to have a negative influence on the time-averaged pressure drop. 

Smith and Swift have experimentally studied oscillatory flow through a nozzle with 
constant diameter, simulating one end of a jet pump.^ In their work, a nozzle is connected 
to open space, establishing a non-conhned jet. A parametric study on the time-averaged 
pressure drop and the acoustic power dissipation is performed, identifying some of the di¬ 
mensionless quantities which describe the flow phenomena: the dimensionless stroke length, 
the dimensionless curvature and the acoustic Reynolds number. Furthermore, a Schlieren 
visualization of the flow held is presented. The formation of a vortex pair and a turbulent 
jet is observed. It is concluded that “extensive numerical studies” are required for a further 
understanding of the minor loss phenomena to control streaming. In a separate article,!^ 
Smith and Swift compare the characteristics of a synthetic (oscillatory) jet to a continuous 
jet in the same experimental setup. The self-similar velocity prohles are found to be identical 
but the jet width of the synthetic jets grow more rapidly than the continuous jets. 

Computational studies related to jet pumps mainly include the work of Boluriaan and 
Morris.ll^*^ In two studies, the minor losses due to a single diameter transition under stand¬ 
ing wave conditions are simulated using a two-dimensional computational huid dynamics 
(CFD) model. The standing wave is generated by either applying an oscillatory body force 
(“shaking” the domain) or by using an oscillatory line source inside the domain. Axial pres¬ 
sure and velocity prohles are presented and the ehect of jetting and vortex shedding on the 
how held is described. The time-averaged pressure drop across the transition is found to 
be a factor of three higher than the quasi-steady solution. In a separate study, a jet pump 
geometry is investigated using a similar CFD model.^ In this case, a combination of two 
line sources with a non-rehecting boundary condition on either side is used to generate a 
traveling wave inside the domain. The how held is calculated for a single jet pump geometry 
and wave amplitude. 

The authors contribution to the held is limited to a preliminary study® where the ehect 
of the jet pump taper angle on the time-averaged pressure drop is investigated numerically 
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and compared against the experimental work of Petcnlescn and Wilen ® A clear decrease in 
time-averaged pressnre drop is observed at higher taper angles, which is one of the motiva¬ 
tions for the work presented here. 


In this paper, the oscillatory flow in the vicinity of a jet pnmp is investigated nsing 
a CFD model which is described in Section Using this CFD model, the performance 
of two jet pnmp geometries with different taper angles are studied. Four different flow 
regimes are described (Section III. A) and subsequently linked to the observed jet pump 
performance. The time-averaged pressure drop and acoustic power dissipation are scaled to 
relate the behavior of the two different taper angles and a comparison with the quasi-steady 
approximation is made (Section III. C). 


II. MODELING 


A two-dimensional axisymmetric CFD model is developed using the commercial software 
package ANSYS CFX version 14.5,^ which has been used successfully in the simulation 
of various (thermo)acoustic applications. *^^*^^ ^ The jet pump is placed in an outer tube to 
study the influence of the jet pump geometry on the flow held. Boundary conditions are 
applied to simulate a traveling wave inside the computational domain; these are discussed 
in Sections II. B and II. C[ In all cases, air at a mean temperature of Tq = 300 K and a mean 
pressure of po = 1 atm is used as the working huid. Three different driving frequencies are 
investigated: 50 Hz, 100 Hz and 200 Hz. 


A. Geometry 


The jet pump geometry is shown in Fig. [T] and is dehned using a reduced number of pa¬ 
rameters: the radius of the big exit Rb, the effective radius of the small exit (the jet pump 
“waist”) Rs,eff, the taper half-angle a and the radius of curvature at the small exit of the 
jet pump Re- Based on these parameters, the other parameters can be calculated. The total 
jet pump length Ljp is 


_ Rb — Rs 

JP — -•) 

tana 

where Rg is the small radius of the jet pump without any curvature applied. 


(5) 
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sin a -|- 1 


cos a 



( 6 ) 


In addition to the jet pump region, the computational domain consists of a section of 
the outer tube on both sides of the jet pump with a radius of Rq = 30 mm and a length 
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TABLE I; Jet pump length Ljp for applied taper angles a. 

(X Ljp _ 

7° 70.5 mm 

15° 35.5 mm 

TABLE II: Dimensions of simulated jet pump geometries. 

Rq 30 mm 

R}j 15 mm 

Rs,eff 7 mm 
Rc 5 mm 


of Lq = 500 mm each for the cases where / = 100 Hz. The influence of the length of this 
section on the jet pump performance and vortex propagation characteristics has been verified 
by comparing with simulations using Lq = 100 mm and no significant difference was observed. 
Although a shorter outer tube section will lead to a reduced computational time, the longer 
length is used in order to study the resulting flow field on both sides of the jet pump in 
detail. For the other two driving frequencies (50Hz and 200 Hz), Lq is scaled relative to the 
acoustic wavelength. 

Two different taper angles, 7° and 15°, are analyzed by changing the jet pump length. 
The corresponding jet pump lengths are shown in Table All the other geometrical parame¬ 
ters remain the same and are listed in Table [H} The dimensionless curvature is Rc/ (27?^ e#) = 
0.36 for both geometries, which is well above the limit for a “smooth” contraction. Hence, 
according to steady flow literature,*^ Kcon = 0.04. The ratio between the small and big 
cross-sectional area is R^^^jj/Rl = 0.22 for both geometries. Because the cross-sectional area 
of both jet pump openings is kept constant, one would expect an identical pressure drop and 
acoustic power dissipation based on the quasi-steady approximation (Eq. [^. Moreover, the 
term (Ag/AbY is small such that the minor losses due to the small opening of the jet pump are 
expected to predominantly determine the time-averaged pressure drop and acoustic power 
dissipation. 

B. Numerical setup 

Within the described computational domain, the unsteady, fully compressible Navier-Stokes 
equations are solved. The ideal gas law is used as an equation of state whereas the energy 
transport is described using the total energy equation including viscous work terms.*^^ No 
additional turbulence modeling is applied as all presented results fall within the laminar 
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regime (see Section III. C). The governing equations are discretized in space using a high 
resolution advection scheme and discretized in time using a second order backward Euler 
scheme. Each wave period is discretized using 1000 time-steps which yields a time-step size 
of At = 1 ■ 10“®s for / = 100 Hz. For each simulation case, a total of Np = 10 wave periods 
are simulated. With a typical computational mesh, the total single core computational time 
is about 40 hours on an Intel Core i7 CPU. 

In order to perform a two-dimensional axisymmetric simulation in ANSYS CFX, a 
computational mesh which extends one element in the azimuthal direction is required and 
symmetry boundary conditions are applied on the originating faces normal to the azimuthal 
direction. On the radial boundary of the outer tube (at r = Rq), a slip adiabatic wall 
boundary condition is used while at the walls of the jet pump a no-slip adiabatic wall 
boundary condition is imposed. 

To generate an acoustic wave, a velocity boundary condition is used at x = 0 which 
oscillates in time according to u{t) = ui sin (27r/f) with ui a dehned velocity amplitude. 
On the right boundary of the computational domain, at x = L, a dedicated time-domain 
impedance boundary condition^ is applied with a specihed reflection coefficient of |i?| = 0, 
as described in Section II. C This ensures free propagation of the acoustic wave without any 
additional reflections being introduced into the computational domain. 


C. Time-domain impedance boundary condition 

A specihc time-domain impedance boundary condition has been implemented in ANSYS 
CFX based on the work of Polifke et This approach ultimately dehnes the pressure 

p{t) on the boundary. The applied pressure is based on wave information from inside the 
domain at a previous time-step which is sampled at a distance Ax = 50 mm from the 
boundary. Moreover, an external perturbation can be introduced on the boundary such 
that any complex reflection coefficient can be specihed. However, in this paper only a non- 
rehective situation {\R\ = 0) is considered. The measured rehection coefficient typically 
ranges from 1 % to 2 %. Further details about the exact implementation and validation can 
be found in the work of Van der Poel which has been carried out as part of the current 
research.®' Additionally, the simulated how held in a tube without a jet pump is compared 
with an analytic wave propagation model®' and excellent agreement is obtained. 


D. Data analysis 

From the CFD results, a transient solution held for all how variables is obtained. In order to 
obtain the complex amplitudes (denoted with the subscript i), a point-wise discrete Fourier 
term is calculated for the specihed wave frequency using data from the last hve simulated 
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FIG. 2: Computational mesh, close-up near jet pump. 


wave periods. The jet pump velocity amplitude is calculated using an area-weighted 

average of the velocity amplitude at the local grid points in the smallest opening of the 
jet pump (see Fig. [^. The time-averaged variables (denoted with the subscript 2 ) are 
calculated by averaging the time-series solution over an integer number of wave periods, 
thus eliminating all hrst order effects. Note that for the time-averaged streaming velocity 
held a density-weighted average is applied: 


U2 = 


(pu) 
(P) ’ 


(7) 


where (...) indicates time-averaging. 


E. Computational mesh 


The resolution of the computational mesh is dehned based on the maximum element size in 
various regions of the domain. A maximum element size inside the jet pump region of 1 mm 
is used which is rehned up to a maximum size of 0.5 mm near the jet pump waist as is visible 
in Fig. 1^ Furthermore, the mesh is rehned near the viscous boundary layer in the jet pump 
region such that a minimum of = 10 elements reside within from the no-slip walls. 

Here 6^, = ^y2fi/ojp is the viscous penetration depth which is 6^ = 0.22 mm for / = 100 Hz. 
At a distance of 50 mm away from the jet pump, a transition to a structured mesh is applied 
to allow for uncoupling of the gradients in the x and r directions and consequently for the 
use of a large aspect ratio as the gradients far away from the jet pump are much larger in the 
radial direction than they are in the axial direction. The axial element size grows towards 
the axial boundaries up to a maximum element size of 10 mm which is sufficient for solving 
the acoustic wave propagation. In the radial direction, a maximum element size of 1 mm is 
used throughout the outer tube. For the current geometry, this yields a total mesh size of 
36 236 nodes. 

In order to validate the computational mesh, the results of three different mesh resolu¬ 
tions are compared based on several key outcome quantities: the dimensionless time-averaged 
pressure drop (Ap 2 , Eq. 14), the dimensionless acoustic power dissipation (AE^i Eq. 15) and 


the propagation distance of the vortex street after ten wave periods, ip. The specified ele¬ 
ment sizes for the three different meshes are listed in Table HTTl In Table. IIVI the total number 
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TABLE III; Specified element sizes for three different meshes used for validation purposes. 
For all other results presented, the medium mesh is used. 

max. element size 

Mesh Nei^BL jet pump jet pump outer 

region waist tube 

Coarse 5 2 mm 1 mm 10 mm 

Medium 10 1mm 0.5 mm 10 mm 

Fine 20 0.5 mm 0.25 mm 10 mm 

TABLE IV: Results of mesh validation study for a jet pump geometry having a taper angle 
of a = 7°. The jet pump waist velocity is = 15.3 m/s, the driving frequency is 

/ = 100Hz. 


Mesh 

^nodes 

Ap* 

AE* 

Ip 

Coarse 

18 442 

0.40 

0.58 

0.81 m 

Medium 

36 236 

0.84 

0.86 

0.92 m 

Fine 

84 618 

0.86 

0.97 

0.96 m 


of nodes is shown. They increase by approximately a factor of two between the subsequent 
mesh refinements. 

Table. |IV| shows the results for the three different meshes using an intermediate wave 
amplitude (ui = 0.8 m/s which yields = 15.3 m/s) and a 7° taper angle geometry 

which is representative for the other simulated cases. The driving frequency is set to 100 Hz. 
A clear deviation is visible for all outcome quantities between the coarse mesh and the 
other two meshes while the results between the medium and hne mesh are comparable. 
The dimensionless pressure drop and acoustic power dissipation obtained with the medium 
mesh, show a difference of 2.4% and 11.1% with the fine mesh, respectively. The vortex 
propagation distance deviates 3.5% with respect to the fine mesh. Hence, it was decided to 
use the medium mesh resolution for all future simulations. 

III. RESULTS AND DISCUSSION 

The described computational model has been used to investigate a range of wave amplitudes 
with the two described jet pump geometries (a = 7° and 15°). It will be shown that the jet 
pump performance can be scaled based on the acoustic displacement amplitude with respect 
to the jet pump dimensions. Defining the acoustic displacement amplitude in the jet pump 
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waist under the assumption of a sinusoidal jet pump velocity as 


Gyp 


\ui,jp\ 

2nf ’ 


( 8 ) 


the two Keulegan-Carpenter numbers can be defined based on the jet pump length and waist 
diameter, respectively: 


KCl = 


Gyp 

Ljp ’ 


(9) 


KCd = 


Gyp 

‘^Rs,eff 


( 10 ) 


Following Smith and Swift iFCi) is similar to the dimensionless stroke length Lo/h while 
KCl is one of the suggested additional dimensionless parameters that may affect the results. 
By investigating two jet pumps of different lengths, it will be shown that KC^ is of high 
relevance to scale the jet pump performance properly. 

Different observed flow regimes will be distinguished and the corresponding flow fields 
will be described. Typical axial prohles of pressure, velocity and acoustic power will be 
described and used to dehne the jet pump performance. Finally, the time-averaged pressure 
drop and acoustic power dissipation will be scaled and shown as a function of the two 
Keulegan-Carpenter numbers. In this way, the jet pump performance will be related to the 
observed flow phenomena and includes the influence of the jet pump taper angle. 


A. Flow regimes 


Independent of the jet pump geometry or frequency, four different flow regimes can be 
distinguished. Examples of these flow regimes are shown in Fig. for the 7° taper angle 
jet pump and a driving frequency of / = 100 Hz. The top graph of each hgure shows the 
instantaneous vorticity field at the last simulated time-step t^ax = 0.1s. The centers of the 
propagating vortex rings can be identihed as local maxima in the instantaneous vorticity 
held. The bottom graph of each hgure shows the time-averaged velocity held U 2 - The black 
line denotes the location of zero streaming velocity. Fig. shows the axial velocity over the 
radius inside the jet pump. The diherent lines are separated (p = 7r/2 in time. Each hgure 
represents a diherent how regime which corresponds to the how regimes shown in Fig. 

In all simulated cases, a vortex pair is formed on either side of the jet pump. How¬ 
ever, for low amplitudes the vortex pairs are not shed and merely oscillate locally with the 
acoustic held. This results in a zero time-averaged pressure drop and negligible acoustic 


power dissipation. An example of this how regime is shown in Fig. ^ where on either side 
of the jet pump a small vortex can be observed. The corresponding velocity prohles inside 
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(a) Oscillatory vortex pair on both sides, no jetting observed. KC^ = 0.09, 
KCd = 0.46, Ap* = 0.04. 



0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 

x[m] 


(b) Propagating vortex to right side, oscillating vortex pair on left side. 
KCl = 0.18, KCd = 0.92, Ap^ = 0.40. 



0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 

x[m] 

(d) Left propagating vortex from waist of jet pump, flow separation inside 
the jet pump occurs. KCl = 1-17, KCd = 5.90, Ap 2 = 0.46. 


FIG. 3: Four different flow regimes are distinguished based on the Keulegan-Carpenter 
numbers KCp and KCd using the instantaneous vorticity helds Vu [1/s] at t = t^ax (top) 
and streaming velocity helds U 2 [m/s] (bottom) around the jet pump for the a = 7° geometry, 
/ = 100 Hz. Black line in streaming velocity held indicates transition from positive to 
negative velocity (color online). 
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FIG. 4: Axial velocity u over the radius at four time instances during the last wave period 
{(f = 0°, 90°, 180° and 270°) halfway inside the jet pump {x = Lq+L jp/2). Figures correspond 
to the four different flow regimes distinguished in Fig. Arrows indicate acceleration and 
deceleration: ► for du/dt > 0 and ◄ for dujdt <0. r/i? = 1 denotes the jet pump wall and 
rjR = 0 is the centerline location. The dashed line indicates the thickness of the viscous 
boundary layer (6^) at the jet pump wall (color online). 


the jet pump are shown in Fig. the prohles are identical but opposite during the forward 
and backward flow direction representing a pure harmonic oscillation. The influence of the 
viscous boundary layer is visible but further away from the boundary a constant velocity is 
observed. 

If the displacement amplitude is larger than the radius of one of the jet pump openings, 
the vortex pair on the corresponding side is shed and propagation starts. Hence, for KCd > 
0.5 vortex propagation on the right side of the jet pump can be observed. An example of this 
flow regime is shown in Fig. 3b where KCd = 0.92. In the streaming velocity held (bottom 
graph), a steady jet in the positive x-direction can be observed. The occurrence of vortex 
shedding corresponds well to an increase in the time-averaged pressure drop. The velocity 
prohles inside the jet pump are shown in Fig. Tb and are comparable to the hrst how regime 
where a harmonic oscillation is observed. 


Fig. 3c shows an intermediate how regime where vortex propagation to the left side 
of the jet pump can be observed in addition to the right-sided propagation. However, the 
how held is still rather asymmetric on both sides of the jet pump which results in a high 
time-averaged pressure drop. The vortex propagation speed is strongly dependent on 
the velocity amplitude in the jet pump waist. Comparing this how regime to the previous, 
the wave amplitude and correspondingly the vortex propagation speed has increased. This 
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results in vortices clearly separated from each other. In the streaming velocity held, a 
recirculation zone inside the jet pump can be observed which is caused by a difference in the 
velocity prohle during the accelerating and decelerating phase. Fig. [4c| shows the velocity 
prohle inside the jet pump at four different time instances. When the huid is accelerating 
{du/dt > 0, indicated by ►), a region near the jet pump wall can be observed where the 
velocity is lower compared to the bulk velocity, regardless of the direction of the bulk how. 
This is initiated when the negative bulk velocity is at its maximum and the acceleration 
changes sign. During the remainder of the acceleration phase, the velocity near the jet 
pump wall (but outside the viscous boundary layer) “lags” the bulk how. When the huid 
starts decelerating (indicated by ◄), a pure acoustic velocity prohle is again observed. This 
diherence in velocity prohles leads to a time-averaged recirculation inside the jet pump. 

The fourth how regime is observed when the displacement amplitude is larger than the 
jet pump length {KCl > !)• An example of this case is shown in Fig. 3d Vortices are 


now displaced from the right jet pump tip through the jet pump to the left, resulting in 
an additional steady jet in the negative x-direction. During the other half of the acoustic 
period, vortices shed from the left jet pump tip are displaced through the jet pump to the 
right and propagate in the positive x-direction, contributing to the existing steady jet on 
this side of the jet pump. These vortex rings are smaller and propagate with a lower speed 
because they are shed at a location where the velocity amplitude is lower than in the jet 
pump waist. After some distance (outside the shown region), the smaller rings merge with 
the larger vortex rings originating from the jet pump waist. 

The existing jet through the jet pump causes time-averaged flow separation inside the 
jet pump which is visible in the bottom graph of Fig. |3dl In contrast with the previous flow 
regimes, now a positive streaming velocity exists close to the jet pump wall. Examining the 
instantaneous velocity prohles inside the jet pump in Fig. ^ gives more insight into the how 
separation process. When the how starts accelerating during the backward how phase, the 
local shear stress at the jet pump wall, becomes zero and the how separation is 

initiated (line marked with ► and negative centerline velocity). After the how reversal, the 
how becomes uni-directional again for the remainder of the wave period. This is independent 
of the sign of du/dt. Although the radius of curvature meets the steady how criterion for 
a “smooth” contraction {Rc/Dg = 0.36 > 0.15), it is expected that the curvature plays an 
important role in the how separation process.^ A detailed investigation of this geometric 
parameter is part of future work. 


B. Axial profiles 

Before describing the relation between the jet pump performance and the jet pump waist 
velocity, the results for a typical simulation case are described for the 7° taper angle jet 
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pump geometry with a driving frequency of / = 100 Hz. These results correspond to the 
flow held shown in Fig. ^ where a steady jet to the right side of the jet pump exists and 
vortex propagation to the left side of the jet pump has just started. On the left boundary 
condition, a far held velocity amplitude of Ui = 0.8 m/s is specihed resulting in a velocity 
amplitude in the jet pump waist of = 15.3 m/s. 

Fig. 5a shows the velocity amplitude and the time-averaged velocity prohle along the 
x-axis. The area-averaged velocity amplitude (dashed gray line) shows nearly incompressible 
behavior where the velocity is inversely proportional to the cross-sectional area. The volume 
how rate Ui to the right of the jet pump is constant to within 0.7%. The irregularities in 
the velocity amplitude at the centerline (solid gray line) to the right of the jet pump are 
caused by the vortex shedding. A positive time-averaged centerline velocity (solid black line) 
to the right of the jet pump indicates a steady jet being formed from the jet opening. Note 
that the jet has only propagated a distance of ip = 0.35 m from the jet pump. Furthermore, 
the steady jet is balanced by a mean how in the opposite direction closer to the tube wall 
(dashed black line) to ensure a zero mean mass hux over the cross-section. 

Fig. 5b shows both the pressure amplitude (gray line, left axis) and the time-averaged 
pressure (black line, right axis). The constant pressure amplitude to the right of the jet 
pump indicates that a traveling wave exists in this part of the domain. On the left side of 
the jet pump a standing wave component is present due to rehection of the acoustic wave 
oh the jet pump surface. This is conhrmed by calculating the rehection coefficient. On the 
right side \Rr\ = 1.16 % while on the left side \Rl\ = 49.9 %. The time-averaged pressure in 
Fig. 5b shows a clear drop between left and right of the jet pump which is one of the main 
measures of performance of the jet pump. While one could simply subtract the time-averaged 
pressure p 2 at two locations on either side of the jet pump, this would lead to inconsistent 
results due to the inhuence of vortex shedding on the time-averaged pressure prohle. This 
is visible in Fig. and becomes more dominant at higher wave amplitudes. Because the 
vortex propagation speed is much lower than the speed of sound, time-averaging over a 
multiple of the acoustic period does not remove the contribution of the vortex propagation 
to the pressure held. Alternatively, a spatially average of p 2 (x) on either side of the jet pump 
is calculated with specihc averaging intervals starting at a distance of 2 ■ from the jet 
pump up to the total vortex propagation distance ip. When no vortex street is present on the 
corresponding side of the jet pump, the averaging is carried out up to the axial extremities 
of the domain. The two resulting spatial averages, left and right of the jet pump, are then 
subtracted yielding the time-averaged pressure drop Ap 2 across the jet pump. For the case 
shown in Fig. a total pressure drop of Ap 2 = 28.8 Pa is generated. 


Fig. 5c shows the axial prohle of the acoustic power E 2 . The acoustic power is de¬ 


termined by the pressure and velocity helds and can be calculated either by direct time 
integration of the transient solution or by using the calculated complex wave amplitudes. 
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(a) Velocity amplitude |ui| at r = 0 (gray, solid)(b) Area-averaged pressure amplitude \pi\ (gray, 
and area-averaged over the cross-section (gray,left axis) and time-averaged pressure p 2 (black, 
dashed). Black lines show streaming velocity U 2 right axis), 
at r = 0 (solid) and at r = |i?o (dashed). 



(c) Acoustic power E 2 . 

FIG. 5: Velocity, pressure, and acoustic power along x-axis using a jet pump taper angle 
of a = 7°. Jet pump waist velocity amplitude is = 15.3m/s and / = 100Hz, 

corresponding to the flow fields in Fig. Vertical dashed lines indicate the exits of the jet 
pump. 
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(a) Time-averaged pressure drop Ap 2 - (b) Acoustic power dissipation AE 2 - 

FIG. 6: Time-averaged pressure drop and acoustic power dissipation for two different jet 
pump geometries as a function of jet pump waist velocity amplitude, a = 7° (•) and a = 
15° (■). Three different frequencies included: 50 Hz, 100 Hz and 200 Hz. 


In the hrst method all resolved higher order effects will be included while in the latter only 
the hrst order acoustics is taken into account. While the integration method will provide 
the most complete solution, the method based on hrst order variables is most relevant for 
thermoacoustic applications as typically only the power transported by the hrst harmonic 
will contribute to a device’s efficiency and any conversion to higher harmonics are considered 
losses.Using the calculated amplitudes of pressure pi and volume how rate t/i, the acoustic 
power is dehned as 

E2{x) = ^iii\pi{x)Ui{x)], ( 11 ) 

where pi{x) is the complex conjugate of the pressure amplitude area-averaged over the local 
cross section. The acoustic power dissipation across the jet pump AE 2 is determined in a 
similar manner as the time-averaged pressure drop. For the current case this yields AE 2 = 
58.8 mW. 

C. Jet pump performance 

By varying the far held velocity amplitude, the relation between the jet pump performance 
and waist velocity is studied. The jet pump waist velocity amplitude ranges from 0.5 m/s 
to 52 m/s which corresponds to acoustic displacement amplitudes from 0.4 mm to 78 mm. 
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TABLE V: Acoustic Reynolds numbers in the jet pump waist and critical Reynolds numbers 
for the simulated cases. 


a 

/ 

max(Re) 

Rcc 

7° 

50 Hz 

486 

524 

15° 

50 Hz 

496 

524 

7° 

100 Hz 

542 

551 

15° 

100 Hz 

521 

551 

7° 

200 Hz 

513 

579 

15° 

200 Hz 

526 

579 


From the velocity amplitude, the acoustic Reynolds number is dehned as 

\ui\6^po 


Re = 


ho 


( 12 ) 


Ohmi and Iguchi derived a critical Reynolds number for the transition to turbulence in 
oscillating pipe flow which can be rewritten to®^ 


Rer = 305 ( — 

0 ,, 


D\^ 


(13) 


The maximum acoustic Reynolds number in the computational domain occurs at the location 
of the jet pump waist. For the simulated cases presented in the following, this maximum value 
falls below the critical Reynolds number which verihes the laminar assumption. Table [V| 
shows an overview of the investigated Reynolds numbers and the corresponding critical 
values. 

The time-averaged pressure drop and acoustic power dissipation for the two different jet 
pump taper angles and three different frequencies are shown in Fig. The acoustic power 
dissipation is similar for both geometries and increases approximately with the cube of 
as was predicted based on the quasi-steady model (Eq. |^. However, for the time-averaged 
pressure drop a large difference can be observed which is a consequence of the change in taper 
angle or frequency. Where the quasi-steady model predicts a quadratic relation between Ap 2 
and \uijp\ (Eq. [^, this is clearly not the case over the entire range of velocity amplitude 
for the simulations. For low amplitudes, a nearly quadratic increase is observed for both 
geometries. However, for the 15° taper angle the pressure drop stagnates and eventually 
decreases at higher velocity amplitudes. Also, the 7° taper angle shows a deviation from the 
theoretical pressure drop prohle and stagnates, but this happens at much higher velocity 
amplitudes compared to the 15° geometry. 
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FIG. 7: Dimensionless pressure drop Ap^ as a function of KCd for two different jet pump 
geometries: a = 7° (•) and a = 15° (■) at 50 Hz, 100 Hz and 200 Hz. Dashed line indicates 
quasi-steady approximation. 


The results at 100 Hz have been compared against preliminary experimental results that 
have been achieved using identical jet pump geometries in a traveling wave experimental 
setup at a driving frequency of 113 Hz. The measured time-averaged pressure drop shows 
the same behavior as the simulation results for both jet pump geometries. The maximum 
deviation is less than 20 % in the region where the measured time-averaged pressure drop is 
signihcant (for Ap 2 > 10 Pa). The experimental results are part of future research and will 
be presented in another publication. 

The effect of the jet pump geometry and frequency on the time-averaged pressure drop 
can be explained by scaling the velocity amplitude using the Keulegan-Carpenter numbers 
based on either the jet pump length or the jet pump waist diameter (Eq. and Eq. 10, re¬ 
spectively). Moreover, the pressure drop and acoustic power dissipation are scaled according 
to-J^ 

8Ap2 


Apl = 


ae; = 


Po\Ul,Jp\‘^ 

?t'KAE2 




(14) 


(15) 


where Ap 2 would represent the difference in minor loss coefficients between the two ffow 
directions assuming the quasi-steady theory (Eq. ^ and would represent the summation 
of the minor loss coefficients assuming Eq. to be valid. 

The onset of the jet pump working is found to be determined by the Keulegan-Carpenter 
number based on the jet pump waist diameter, KCd, and is shown in Fig. for the two 
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(a) Dimensionless pressure drop (b) Dimensionless acoustic power dissipation 

AEI. 

FIG. 8: Dimensionless pressure drop and acoustic power dissipation as a function of KCp 
for two different jet pump geometries: a = 7° (•) and a = 15° (■) at 50 Hz, 100 Hz and 
200 Hz. Dashed lines indicates quasi-steady approximation. 


investigated taper angles and three different driving frequencies. No time-averaged pressure 
drop is measured for low values {KCd < 0.5). This corresponds to the flow regime described 
in Fig. 1^ where there are no minor loss inducing flow phenomena observed. In the low 
amplitude regime, the effect of the jet pump taper angle is negligible and both geometries 
follow the same line. However, at higher amplitudes the curves become deviant and the effect 
of jet pump taper angle becomes apparent which is well accounted for by using KCp rather 
than KCd- Note that the effect of frequency is well accounted for by using the acoustic 
displacement amplitude as a scaling parameter rather than the velocity amplitude. 

Fig. shows the dimensionless quantities Ap 2 and as a function of KCp. Both 
jet pump taper angles show the same trend indicating that the jet pump performance is 
directly related to the jet pump length. A maximum in the dimensionless pressure drop 
is visible around KCp ~ 0.5 for both geometries. As soon as the displacement amplitude 
becomes large with respect to the jet pump length, the total minor loss coefficient begins 
to decrease again. This corresponds well to the observed asymmetry in the flow fields as 
discussed in Section HI. A At KCp ~ 0.5 the flow field is highly asymmetric between both 
sides of the jet pump (Fig,3c) while for KCp > 1 a large part of the asymmetry is lost and 
flow separation occurs. The time-averaged separated flow reduces the “effective” diameter of 
the big opening (Fig. 3d). This will consequently lead to more symmetric minor losses and 
hence to a decrease in the time-averaged pressure drop while the acoustic power dissipation 
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will still increase. In this regime, the quasi-steady approximation is not valid anymore and 
an adjustment is required. 

Comparing the obtained maximum value of Ap 2 with the quasi-steady approximation 
(dashed line) in Fig. 8a, a close match is observed for both taper angles. This suggests the 
quasi-steady approximation can be used to predict the time-averaged pressure drop but only 
for the optimal situation where KCl ~ 0.5. The acoustic power dissipation is predicted well 
for all cases by the quasi-steady model as shown in Fig. 


8b 


The assumed cubic relation 
between AE 2 and |'Ui,jp| is conhrmed by the fairly constant value of AE 2 , especially for 
0.5 < KCl < 1-5. In this region, the jet pump acts as a pure acoustic resistance and its 
behavior is not affected by a change in flow regime. 


IV. CONCLUSIONS 

A computational fluid dynamics model is successfully used to simulate the oscillatory flow 
through two different jet pump geometries under traveling wave conditions. The relation 
between the obtained time-averaged pressure drop and the acoustic displacement amplitude 
is investigated together with the acoustic power dissipation. 

Four different flow regimes are distinguished based on the Keulegan-Carpenter numbers 
and the observed flow phenomena are related to the jet pump performance. For KCd < 0.5 
no vortex shedding is observed, resulting in a negligible pressure drop. At KCl ~ 0.5, the 
highly asymmetric flow held leads to a substantial pressure drop and a steady jet to the 
right side of the jet pump is observed. The measured jet pump performance in this regime 
corresponds well with the quasi-steady approximation. 

When KCl > 1, vortices are shed through the jet pump to the left resulting in an 
additional jet to the left side of the jet pump and how separation in the jet pump is observed. 
Reducing the asymmetry of the how held consequently leads to a decay in the time-averaged 
pressure drop but increases the acoustic power dissipation. In this how regime, an adjustment 
to the quasi-steady approximation is required to reliably predict the jet pump performance. 

Several additional geometric parameters have not been considered in this paper and 
additional research is required in order to distinguish more precisely between the diherent 
how phenomena. This should provide insight into whether the decay in Ap 2 for high values 
of KCl is solely caused by the jet pump length or by the jet pump taper angle as was 
postulated in previous work.® Moreover, the inhuence of the radius of curvature on the jet 
pump performance is not investigated here. The radius of curvature is expected to have an 
ehect on the how separation inside the jet pump and consequently on the overall performance. 

A hrst step towards a better understanding of the physics behind jet pumps is made 
but more attention is yet required to reliably predict a jet pump’s performance. A thorough 
numerical parameter study on the various (two-dimensional) geometric parameters together 
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with experimental research on three-dimensional geometry variations are the next steps to 
be undertaken to provide insight into to the jet pump scaling problem. 
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